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ABSTRACT 


The Haltiner theory of saturated steady state 
CGmvectivescummenecewhicheineiudes the transfer of heat 
and momentum by lateral diffusion as well as the 
entrainment of environmental air is discussed. A 
modified model is presented which includes, as a 
separate term, the drag of the liquid water carried 
along by the ascending current. The system of equations 
is integrated numerically for several cases previously 
considered by Haltiner allowing entrainment to assume 
negative values (detrainment) and the results are 
compared. The results are in good agreement, and show 
tiene meGeetfect of the Opposing actions of liquid 
water drag and detrainment is a reduction in the 
SoMoleeo Veruvcals veloc Fics. This reduction appears 
Opes mee ch noms inm tne righty direction -on the 


basis of a comparison to some observed data. 


[nesWwerverswiches TO express his appreciationwrern 
the assistance and encouragement given him by 
Professor G.J. Haltiner of the U.S. Naval Postgraduate 


school in this investigation. 
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1. Mite OCucéon 

The concept of simple "adiabatic parcel ascent" as 
eat irst Approx rmavtionelowcen eet ivememenomena has: long 
Been sin usess itehas, nowever,.aaumeer ot shortcomings. 
As a result of evidence (Stommel,|s] ) that a rising 
@anvectCive current acts like a jet, dragging in surround- 
ing air as it progresses upward, Austin and Fleisher 
Suggested that this process, called "dynamic entrain- 
ment", is necessary to satisfy mass continuity. (1] 
Haltiner{4] summarized the work of Houghton and Cramer,lL5] 
Bunker,{2] and Malkus,|[¢] and suggested a theory of 
saturated convective currents which includes the effects 
of lateral diffusion of heat and momentum, and 
extended computations to elevations typical of large 


cumulus clouds. 


Haltiner [4] derived a system of equations which 
eescrive, a savuravecmoomveccive current. The derivagien 
of these equations is essentially the same as the 
derivation of the modified system of equations presented 
below. The Haltiner model assumes that during 
detrainment the cross-sectional area of the rising 
column remains constant. However, when mass continuity 
considerations calls for the ejection of mass 
(detrainment), a second system of equations is applied 
based on the assumptions that the rate of entrainment 


remains identically zero and mass continuity is 


N 





satisitile@ by an increasime crosse-seetional area in 

the rising air. In this paper the detrainment phase 
will be treated by retaining the assumption that the 
cross-sectional area of the rising air remains constant, 
thus permitting the entrainment to assume negative 
values. 

The effect of the drag of the liquid water carried 
along in the rising air was not explicitly considered 
fa thesiaiiremersmodel, aitviemwem it can be qualitarively 
iclided@iye ad UuStine = wte rave Of diffusion of 
momentum. Here the drag effect of the liquid water will 
be specifically included as a separate term. 

The following notation will be used throughout 


this Paper : 


a Tempe aac ure 

Z height 

q Sveeil 1c Numidasy 

Cp specific heat of air at constant pressure 
L latent heat of vaporization 

g gravity 

Ww velocity 

k emi tusvon "eeerfictent 

ei gas constant for dry air 

Ry gas constant for water vapor 

N atl 

M Mass rate of flow per unit time 





eae 


"specific humidity" of liquid water 


cross-sectional area of rising column 
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density 
p pressure 
¥ lapse rate 
subscripts: 
e enviroument 
V virtual 
O init@val value 
Ie, The system of equations 

ao. The continuity equation. 
Consider a vertical current of air in which the velocity, 
Wy is horizontally uniform. If the field of vertical 
velocity is not uniform w may be considered to be the 
mean value in the horizontal at any given level. The 
mass of air per unit time flowing past any given level, z, 
is given by M = ow G 
where M is the mass, -~ the density and 0 the cross- 
sectional area of the vertical current. Under steady 
state conditions the mass at an adjacent level # +dzis 
the sum of the mass rate of flow at level = , and air 
entrainment fromz to z+dz. Thus, 
(1) aM do Spain 2 


mM ia w on 


wrals & '«¢ > = 


re? ee ee “a” ‘=> 

a “=e 2 FF eu 
ee (= © eee 
eel Eel <<) ae + = 





Defining a parameter N= 1inM, we may write 


c|N 

dz amelie w cl2 G clz 

Dif ferentiation of Chexequationsoimstate p= ty 
yields 


¢ 


Jl 


We next make the assumption that at any level the 


pressure in the rising air is the same as that of the 


environment, then dp — BE Dit 
P dz P az 

From the hydrostatic equation @p/foz= —- 42 abt 
C2 


follows that 


a SP L Vrele es: 
P dz Ne OZ ays 


We may now eliminate the density term from (2) by 
substitution from (3) and (4) to yield the continuity 


equation for the rising air 


(eed | Giedvtetkdie., Soe edd 


a 





dz wdz g dz Rt. 80 dz 

Cee Ne weauatlomno: Motion. 
The changing mass of the rising air due to entrainment 
or detrainment requires the use of an equation of 
motion for a variable mass. If a mass M is moving at 


velocity w at time t combines with a second mass aM 


moving at velocity w, during a time interval dt 


Ly 





and w+dw represents the velocity of the combined masses 
av Uaioe + Oty CHEN Gueemenenunmoma. Lime t and t +d 


are related byeimercquation 
(MrdM\(Wwrdw) = Mw rte Mdt tw adM 


Here E represents the external forces per unit mass. 
This relationship yields the equation of motion 
eee ee yt ah 
dt Misdit M dt 
ne extemmaletorces acting on the system inclimemrne 
pressure and gravity forces, which may be expressed as 
g(T- Tye) / Te , and frictional forces, 
Enewtattermenciuding the drag of the liquid water@egmried 
by the ascending column of air. The equation of motion 


for the liquid water rising with the air may be written 


Ciiucdw a flac 
ge B- Ag +d 


where | is the mass of liquid water per unit mass of 

air which we will term the "specific humidity" of the 
neaveramamManWemepresents the drag of the air acting en 
ene liquid water of mass // mm Ine negative of D isa. 
drag of the liquid water on a unit mass of air. From 


eq (7), we have 


—|) = —}q - dy 


substituting in eq (6) we obtain 


dw = -~-dw He ae cit ae nes W - WW ad M 
ze = ge - 4a +g (Se) — wowed. F 


Simplifying we have 
; ee el ae i . 
Wal) dw 2 g (hole 2) - (w-weidN 4 OF 
dt dt 


Me 
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assuming steady state conditions we obtain 


- : wT, ; 
(a) Wa = | 9 toe 8) iwi rll 4 F 


dz (1+) We dz 





c. Lateral diffusion of heat and momentum 

In addition tosthe exechanse due vOvemirainmenct, 
turbulent exchange of heat and momentum must also be 
considered. Priestley €7J(1956) investigated the 
mechanisms of diffusive heat and momentum exchange in 
a plume of hot gases. Haltiner [4] suggested that, 
excepl forezne scale of the phenomena, the diffusive 
Exmclanmpe preccsses in a convective current are Ss iimelar 
to the processes defined by Priestley. The effects of 
Paeral aii iusiton were included in the Haltaner modem 
Pemepne extenvion of the Priestley equations to conmvecuive 


currents. Yhem™neat loss per unit time is taken to be 
(9) -~hk, o,f (@-@ 
C 
or approximately 
mOv—=— ay C(t |, ) 


Itv-as also necessary to saturate the environmental’ air 


entrawmiecd Ante the updratpmand a” term of the form 
(11) ~k, L (4- Ce) 


iomGcgumcdgm@unere Lois the latent-heat’ of Vaperization. 


The loss of momentum per unit time per unit mass is 


6 





Seellar’ vo thay Ge vio 1Oca er cae and is assumed to 
be proportlonad bate the dei ierceicemeewween the velociry 
of the convective current and the environment, with the 
latter taken to be zero, with the result ~k,w where a, 
has the dimension second™, Wartho this formmiorethe 
friction term, eq (8) becomes 


a2) BY = din [S(TRA)- mand — he] 


CA (14 4) fe ve. ra Fae 


the dimensions of &, imply that diffusive losses are 
imd@éependéent™ of updraft velocity. ™ Haltvimer “sugeecrs a 
second approach to the problem of diffusion which makes 
the loss of heat and momentum dependent upon updraft 
vele@city. The first approach, showever, will be employed 


in this paper. 


d. Thermodynamic equation. 

In describing the thermodynamic process undergone 
by the rising airya number of factors must be incitdeds 
The following must be taken into account in accordance 
With the First Law of Thermodynamics: the heat released 
by condensation, the heat required to saturate the 
Elubearned environmental air and to Ralse 10s (temper ure 
bor Uniav Of thesrising column, the Nest est Tomrnec 
environment through diffusion, the internal energy 
changes and the work done during expansion. The 
Besulting equation embodying these iectors is as 


follows: 








ice ee SI ~i(q-qe)aM ~ <> (T-TlM 
— Mk, fe,(T-%) +L (4-qe) Jat = 
M Bec ae RT, d(inp) | 

The terms on the left, in the order of their 
occurrence, represent respectively, (a) the heat 
released in condensation; (b) the heat required to 
saturate the entrained air, (c) the heat required to 
bring the entrained air to the temperature of the rising 
column, (d) the heat lost through diffusion in the 
interval dt required for the air to rise a distance 
ize (heewienG esteem or the equation, represents tie 
combined internal energy change and work by the gas in 
expanding from pressure p to pressure p+dp. By 
substitution from eq (4), eq (13) becomes, with 


appropriate manipulation, 


—Ldq ~ [c,(T-R)+L(q-4q)|d 
dz 


dz 


—di\ oc, (T- Te) + & (q-4e) 
| aan : 
= ct gf 


Ww 


IEF) 


— + 
Bice Wee 


e,. Moisture and liquid water. 





it has been assumed that the rising air is sat-= 
urated; the specific humidity of the rising air is then 
a Tunctionmel temperature and pressure. As a good 
approximation we have the expression gq = 6219p where © 
is the saturation vapor pressure. By differentiation 


this becomes 





wer —e eee 


(15); dq _ 14 de i dp 


— 


Cf dz Gigcle ia AZ 


The last term in eq (15) may be replaced with 97K, Ty. 
pee 





from eq (9), and from Clapeyron's equation, = ie. 
where RL is the gas constant for water vapor, an : 
alternate expression for the second term is obtained, 
yielding 

& 
Roleese 


cr 





(16) 


Bay 


He eased 
gis 


el 


Ther ligu@de water content of the cloud is the amounteol 
moisture condensed less the amount reguired to saturate 
pac envir~ommenoal Aly Entrained and that exchanged 


through diffusion. Defining ¥ as the "specific humidity" 


of the liquid water we may write 
(17) dim) =- Mdag ~(q- Ge) dM - MR, (4-de) dt 


Expanding eq (17) yields 


aadt, 42 4 (+q-q) dn 
ale 


dz dz Sec) 0 


Equation (17), however, is not valid during detrainment 
and the change in liquid water content must be piece- 


Wise defined. For the detrainment phase we have 


(19) Jd(M2) = —~ Mdd +ldm- Ma, Cq-Ge) At 


whick upon expansion leads to 





(20) GX dy ky (gq -4e) 
dz Az W 
Two additional relationships between temperature and 


Veruual bemperature will be required. sone of these is 


the well known relationship 
eee + .61q) 7 


the second relationship is obtained from eq (71) by 


@aeiferentiating with respect to Z; this yields 


a (it -6IQ\dT + cit dg 
clz Jz az 


Substituting from eq (16) this may be written 


ee a (alba »61gG + ight \dt Tag 
dz - eae cl Z ee 


By using some mean values this may be approximated by 





-) a (| + 12-5q)ot 
aa 1) dz 

{, Tne environment. 

The environmental parameters |,, ly. , de, andw, 
appear in one or more of the equations derived above; 
it is therefore necessary to define the environment in 
some way. We have wide latitude in choosing the struc- 
ture of the environment; 1% may, 10r instance, scommerm 
to an actual sounding. A simple model was used in this 


investigation. The lapse rate %& and the relative 


humidity were assumed constant, and W, was assumed to 


10 





be zero. Under these assumptions Ts may be expressed 


as 
(23) [| = We ye 


dime environmental specific humidity may be determined 
as follows: From eq (16) the saturation specific 


MPumidity 1s given by 
OC go = @ L Fsx. %e a 9 dsc. aez 
i oe sae 


and since the relative humidity is constant G.=Fq., we 





may write 


de. (~ ater. + Bie jaz 


(24) RO ie re 


Me the lapse rate ca Clas bigs abel CemioemecurnGe ca. » and 


the relative humidity are chosen, then the values ofl,, 


Tye 


, and q, may be obtained at any level from equations 


Cem, (23), and (24). 


BA The complete model. 

Having defined the environment in sub sotticam 
Pemations (5), (12), (14), (16), (20) ane (2) ora 
constitute a system of six equations in 7 unknowns,w, 
iaae 1qads N, andg . We must now establish a 
melationshbip defining ¢ . Ihe stiuplest solutionmeana 
the one which will be employed here, is to assume that 

0 is identically zero. This assumption is quite 
ively an Ooversimplification,. bul on the otnNer hand 


presently available information on the behavior of the 


ea 





cross-sectional area of an updraft is scanty and somewhat 


com: licting: 
The system of equations (5), (12), (14), (16), (21), 


and (18) or (20) may be solved first for the cloud 


temperature lapse 


= ay (ype egret." (Tye eA) ~ 3 Bs - de | 


+ by | (4 ~Ye) 4 Ce -_ ak + 49 
w c ye | Rate 


gare (alias) (=a) + $207 TE) ] Y 
. ee L 


\ 
——— aa : 
i le 2a ly 

and then successively for the other parameters 


dN Pee Ey NT Wve mi Y i od ; 
(26) O27 an aki t| tu — oop - +30] 


ve dive 


Mle 





dl : i -~“G dN _ h a (as ee teu 
(27) d on | Ts 7 “\ Bice 9] 
d L d I 
(28) 0 qg 


| = Aout 
dz RT? dz Bele 


(29a) IA _4q _ (4+49-4.) aN = &, (4 - qe) 


dz dz dz Ww 
or 
29) dk dq _ ey (4q- qe) 

dz cl Ww 


Equations (21) through (24) are used as needed. Given 


initial values of T, N, w, 9 , »g» g., the foregoing 
system of equations may be integrated by numerical 


towhneds over the desirederange of Z. 


i 





oerce™adequate "data concermime the vertical welocity 
of the entrained environmental air is lacking; eq (27) 
will be further simplified By assumime that w,=0. 
Moreover since typical values of are on the order of 
.005 we may place 1f1+4) + 1 . The error introduced 
by the latter approximation will be on the order of $ 


of one percent. 


If eq (27) is modified as indicated above, it becomes 


ale "I A) 


= 


(30) 


| 
| 
' 





Ceien 
mM 1S 
< 
Cie. 
m !8 
sx 
< fo 


a (orm which differs from the corresponding equalmaenm 
in the Haltiner model only in the last term, which has 


the form 


a lve 
Similarly in the Haltiner equations for di/dz anda 


a 


dN/dz, ee incre a (%e — Ne is replaced by a 
We Flee Ww We 





in the modified model. Eguations (28) and (29a) are 
exactly the same as their counterparts in the Haltiner 
model. The only other difference in the two models is 


the appearance of the additional eq (29b). 


Ly Computations and Results. 

The foregoing system of equations was integrated 
by numerical methods on an electronic digital Computer, 
the National Cash Register 102A. With a given set of 


initial values, properly scaled inputs to the integration 


Le 





meweine were computed for the model given above. These 
Valuesswere theneintegraved over -anm@interval of -leemeters 
yielding a new set of “initial" values for the level z+ 
10 meters which served as an input for the computation 

ot a néw sebeot inputs for ne intesrariom routine. 
Computation proceeded in this stepwise fashion until the 
top of the cloud was reached. 

Computations were made for four sets of initial 
eonditions. in eacn case Ene snitial values Ceommesponaed 
fo one of the sets of initial values used by Haltiner|4 ] ° 
In Case (a) and Case (b) the same set of initial values 
were used, but in Case (a) only the effect of detrainment 
tee COnsiagercd, “Ihe eirect =o: the drap of the Itaberd 
water was eliminated from Case (a) by an appropriate 
modification to the complete program. In Cases (b), 

(c), and (da), both the effects of detrainment and the 
drag of the liquid water are included. 

The initial values used in the computation are as 
ol Lows : 

Case (a) 
Ty = 20C, Teo = 19°C, wo 71 meter/second 
ka =__.001/second, Ke =.7  C/100 meters 
r = 80% 
Case (b) 
Same as Case (a) except that both detrainment and 


water drag are included. 
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Case (c) 
Same as Case (b) except that = .6° C/100 meters 
Case (ad) 


Same as Case (b) except that r= 90% 


The results of the computations are presented in 
figures 1 through 6 and table l. 

Figures 1 and 2 present the results of Cases (a) 
and (bv), in which the parameters w, 4 ; M/Mo and A T 
are plotted as functions of height. The parameter M/Mo 
is the cumulative change im mass experienced by amgunie 
of air due to entrainment. A value of two, for example, 
midicates that the mass hes doubled. ihe, paramever= 7. [ 
Poecct Medeas 7] = hae. 

Figures 3, 4, 5, and 6 present a comparison of the 
rpesults Obtained by Haltiner|"4]for the set of values 
listed as Case (a) with the results of Cases (a) and (bd). 
Each parameter is plotted separately as a function of 
heiehnt. The values computed by Haltiner and the values 
from Case (a) should be the same at and below the level 
Wiegemcdcunainment beginsa Théessystematic difference aa 
two values below this point is the result of minor 
differences in the computational procedures (primarily 
in the computation of the environmental parameters) 
followed in this investigation and those employed by 


Lemay Inetweweeinis systematic iWitterence 1s noteewanaaacant 


an 





mmc self, but itedoees illustrate how small differences 
at one level propagate upward and tend to accumulate. 
Tne el lects ol derrcaimmnentwanimene crag of Vite 
liquid water are best illustrated in figures 4 and 5. 
Tt is apparent that the effect of detraimient alone is 
itemeaccelerate the verticaimcurrent.. [ne choudmnei ent 
is greater in Case (a) than it is in either Case (b) or 
the corresponding case comcuted by Haltiner\|4]. It is 
also apparent from the figures that in Case (a) the 
magnitude of the maximum values of w and M/M,and the 
elevation at which these maxima occur are greater than 
emener the values obtained by Baltime. or the vatnes 
of Case (bv). These results, and the results of Cases 


(c) and (d) are summarized in table one. 


os Conclusions 

The number of cases investigated was not exhaustive, 
buy Wasmem ticient to show that tre moditied model 
yLelds results that conform quite closely to the resuilgs 
of Haltiner| 4]. In particular, the cloud heights were 
substantially the same but the vertical velocity and the 
amount of entrainment are reduced. These results are 
the natural consequence of the combination of detrainment 
and the drag of the liquid water. Detrainment on the 
one hand provides an acceleration, while the drag of the 
itewia Waver acts as a retarding orce ]  inesem e.ces 


Pema vor balance one another, but the drag verethe liquid 


16 





Beaver Somewhat overbalances the acceleration due to 
detrainment. 

Haltiner L4J] investigated an actual cloud for which 
observational data was available, and found that his 
results were in good agreement with the obDserved values. 
He did find, however, that his values of vertical 
velocity, although of the right order of magnitude, were 
somewhat in excess of the observed values. To the 
extent that the modified model reduces vertical velocity, 
mac mOdil ications appearm@ee be a2 Ceomeect ion in themriene 
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